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ABSTRACT 


The statistical properties of the radar derived rainfall obtained 
during the GARP Atlantic Tropical Experiment (GATE) are used to 
derive quantitative estimates of the spatial and temporal 
sampling errors associated with estimating rainfall from bright- 
ness temperature measurements such as would be obtained from a 
satellite-borne microwave radiometer employing a practical size 
antenna aperature. The analytic results of this study provide a 
basis for a method of correcting the so-called beam-filling 
problem: i.e., for the effect of non-uniformity of rainfall over 
the radiometer beamwidth. The method presented employs the 
statistical properties of the observations themselves without 
need for physical assumptions beyond those associated with the 
radiative transfer model. The simulation results presented offer 
a validation of the estimated accuracy that can be achieved and 
the graphs included permit evaluation of the effect of the 
antenna resolution on both the temporal and spatial sampling 
errors. For example, it is found that the bias associated with 
the antenna having a resolution of 32 km can be expected to be 
corrected to within ±6 percent of the monthly mean at a one sigma 
confidence level with two observations per day over a 256km x 
256km region for which the GATE data can be considered typical. 
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X.O INTRODUCTION 


The Electronically Scanned Microwave Radiometer (ESMR) flown 
on the Nimbus-5 spacecraft demonstrated a capability to sense 
rainfall over ocean areas. The purpose of this study is to 
provide an increased understanding of the potential capabilities 
of microwave radiometry for remote observations of rainfall and 
to develop quantitative estimates and corresponding confidence 
bounds of the measurement accuracies that can be achieved. 
Results presented here should be of value in establishing the 
range of scientific requirements for remote observations of rain- 
fall that can be met with current technology and «:hey should be 
of further value in designing any future flight programs invol- 
ving ESMR-like instruments. 

The primary objectives of this study include: 

1. Determination of the extent of bias errors due to non- 
uniform and incomplete filling of the radiometer beam- 
width (i.e., beam filling errors) resulting from direct 
conversion of observed brightness temperatures to rain- 
fall through radiative transfer model. 

2. Development of the relationship between such bias errors 
and the antenna resolution as a function of brightness 
temperature parameters. 

3. Development of a method based on this relationship and 
estimates of the brightness temperature parameters 
derivable from the observations for unbiased estimates 
of the mean rainfall. 

4. Assessment of the quality of this method in terms of 
expected accuracy and confidence levels, 

5. Computer simulation to validate that the bias correction 
method developed performs as expected. 
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Tlie GATE* radar observations during the summer of 1974 proved 
to be nearly ideal for this study. They were used to derive the 
statistical properties typical of rainfall that were essential to 
this study and they were also used in the computer simulations 
after converting the radar measurements of rainfall into corres- 
ponding brightness temperatures using the radiation transfer 
model of Wilheit, et, al. [1977J. The details concerning con- 
struction of the brightness temerature scenes are presented in 
Section 2. The beam-filling problem is presented in Section 3. 
The formulas and relationships necessary to attain the objectives 
of this study are presented in Section 4. A discussion of the 
simulation results is presented in Section 5 and a summary is 
given in Section 6, 


* GATE (Global Atmospheric Research Program) Atlantic Tropical 
Experiment. 
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2.0 RADIATION TRANSFER MODEL AND PREPARATION OP THE SIMULATED 
DATA 

This section briefly describes the radiation transfer model 
of Wilheit, et, al, [1977] as it applies to the estimation of 
rainfall from brightness temperature observations and to our 
conversion of the GATE radar data into simulated brightness tem- 
perature scenes. 

Rate a nd Brightness Tai^perat ure Relati onship 
The thermal microwave radiation from the ocean surface as 
seen from space is modified by the liquid water and water vapor 
within the intervening atmosphere. The eraissivity of the surface 
depends on its dielectric constant which for ocean surfaces 
varies approximately inversely with water temperature (T^^), 
However, at an incidence angle of approximately 55° and in the 
region of 19 Ghz, these are nearly compensating factors. As a 
result, the brightness temperature (T 3 ) which is the product of 
emissivity and water temperature is nearly constant for smooth 
water. The wind at the water surface can effect the emissivity, 
but this effect is relatively insignificant. Therefore, the ocean 
presents a nearly uniform background radiation source so that 
variation in the received radiation results from variations in 
the attenuation of the radiation that are mostly due to varia- 
tions in the liquid water content of the atmosphere. 

The equation of radiative transfer through the atmosphere 
given by Wilheit, et. al. [1977] is 
dTj,(G) 

— S aT„(e) = S 

dZ ® 


Tg( Gg )f( 0 jGg ) Sin Gg dGg+aT^ (2.1) 
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where? 

Tg(e) » radiance (or brightness temperature) in the direction 
° specified by the polar angle 0, 

a * total attenuation coefficient, 

S « scattering coefficient, 

a « absorption coefficient (o'«a‘fS), 

« thermodynamic temperature of the atmosphere; and 

F(e,0e) ^ scattering distribution function for scattering from 
angle Ogtoe. 

In 2.1, the absorption is attributed to all causes including 
gases, water droplets, water vapor, and ice crystals. At wave- 
length near 1 to 2 cm, only molecular oxygen, water vapor and 
water droplets are taken into account. However, contributions to 
the total absorbtion due to molecular oxygen and water vapor are 
relatively small and assumed values are used. The contribution 
of water droplets was calculated using Mie theory and a raindrop 
size distribution given by the Marshall-Palmer distribution, 

N(r) « Nq e-«»^ (2.2) 

where, N(r) is the number density of droplets per unit size 
interval, r is the radius of the droplets in centimeters, Nq is 
0.16 cm""^, and 

6 = 81.56 R-0.21. (2.3) 

In 2.3, R is the rain rate in mm hr”^. Equations (2.2) and (2.3) 
provide the relationship betv/een N(r) and R, needed to solve 
equation (2.1) for Tb as a function of R. 

The relationship between Tg and R was solved by Wilheit, et. 
al, [1977] for a layered atmosphere assumed to have v;ater drop- 
lets from the freezing level to the surface. They modeled the 
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watec vapor content by the specifying a linear variation in 
humidity from the surface to the freezing level and they assumed 
that cloud water droplets have « density of 25 mg cm*"^ for the 
first 0.5 km below the freezing level. Results of their calcula- 
tions with an assumed lapse rate of 6.5°C per km are displayed in 
Figure 2-1 where the freezing level has been left as a free 
parameter. It should be pointed out that their model indicates 
that the brightness temperature increases monotonically with rain 
rate only up to 20 to 50 mm hr*“^ after which the back-scatter of 
larger rain-drops causes a decrease in brightness temperature. 

2.2 Simulation of Brightness Tr'^perature Data 

The radar measurements of rainfall, collected during the 
summer of 1974, over the GATE B-scale area, located approximately 
1000 km off tlie west coast of Africa, are described in detail by 
Hudlow [19793. The radar data used was specifically prepared for 
this study by the NOAA {National Oceanic and Atmospheric Admin- 
istration) which provided the processed digital data (for Phases 
I and II) for each 4km x 4km rain cell at the original 15 minute 
observation intervals. A description of this specifically pre"- 
pared GATE radar data is given by Laughlin and Gupta [1980]. It 
is recognized that although the radar measurements of rainfall 
may be somev;hat erroneous, this will not have any significant 
effects in this study as long as the simulated ESMR has resolu- 
tion larger than 4 km and has sampling interval larger than 15 
minutes. 


5 




ORIGINAL PAGE IS 
OF POOR QUALITY 



6 


Figure 2-1 

Calculated Brightness Temperatures at 1:55 CM as a Function of 
Rain Rate for Melting Levels of 1, 2., 3, 4, and 5 KM (from Wilheit, et al., 1977) 


Bach pceolpitahion field consisting of an array of lOOxlOO 
cells was mapped into a corresponding brightness temperature 
field using the transformation curve of Figure 2-1. This 
transformation was done for each individual 4km x 4km cell in the 
precipitation (rain) field, A freezing level of 4km was assumed 
without any loss of generality. For computational ease, the best 
analytical fit for this transformation curve was found and has 
the following form^ 

Tg » 271 - 107 exp(-0.182 R) for R ^20 mm/hr (2.4) 

« 271 - 0.1944 (R-20) for R >20 mm/hr 

This is a nonlinear relationship between R and Tg and, as will be 

seen later, introduces some error while estimating mean rainfall 
from observed brightness temperatures. These simulated bright- 
ness temperature fields formed the basis and provided the data 
for all the analysis performed in this study. 
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3.0 ERROR IKi ESTIMATING RAINPALIi PROM BRIGHTNESS TEMPERATURE 

The accuracy with which passive radiometers can estimate 
rainfall from observed or simulated brightness temperatures 
depends on several factors. Tliese include the non-linear and non- 
singular behavior of the Tq vs, R relationship, temporal sampling 
rate, and the spatial resolution of the antenna. The behavior of 
Tq and R relationship reflects the instrumental errors such as 
biasing of rainfall measurements due bo incomplete filling of the 
radiometer beamwidth. These problems are addressed in this 
section. 


3.1 Beam Filling Problem 

Non-uniform and incomplete filling of the radiometer beam- 
width in conjunction with the fact that Tg vs, R relationship is 
rsoi -linear and non-singular introduces a bias in the estimate of 
mean rainfall and this bias causes the rainfall to be consis- 
tently underestimated. 

To understand this problem, let us examine a f ield-of-view 
(fov) of sufficiently large area such as 64km x 64km formed from 
an array of n cells with the cell size being much smaller in area 
than that of the fov. Assuming that the radiometer antenna has a 
resolution equal to cell size, the radiometer v/ill measure the 
average brightness temperature, Tt>i , in cell i for each cell in 
the fov. 

The mean rainfall, over the fov is then obtained accor- 
ding to equation (3,1). 


R. 


_L 

n 


n 

i«l 


(Tbi> 


(3.1) 
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Where, represents the inverse transformation for converting 
brightness temperatures (T^i) into corresponding rain rates {tj), 
fhis is easily obtained from equation (2.4) and expressed as; 

t~? r, » -i- «n / (3.2) 

where, a=271.0, b«107.0, and c«0.182 are constants, 

/ 

In order to determine the existence of bias in estimating 
rainfall from observed brightness temperatures, let us now scan 
the same fov with another radiometer having an antenna of lower 
resolution such as k times the cell size, k being a positive 
integer greater than 1. The mean rainfall, R 2 , over the same fov 
can be obtained according to equation (3,3), 

jk \ 

^ (3.3) 



The existence of bias will be positively confirmed if the quan- 
tity - R 2 is positive. In addition, the value of this quantity 
would indeed be the amount of bias error. To obtain a simple 
expression for R^ - R 2 r let us substitute equation (3,2) in 
equations (3.1) and (3.3) and simplify. 
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Using a series expansion: Rn(t-m) 
and simplifying: 



} 



V7here, represents the empirical moment of brightness 

temperatures observed with antenna having a resolution equal to 
cell size and mp^j^ represents the p^^ empirical moment of bright- 
ness temperatures observed with antenna having a resolution equal 
to k times the cell size, A close examination of equation (3.4) 
reveals that each term in the summation is positive and that the 
entire series is convergent. Furthermore, it is clear tha as k 
increases each terra in the summation also increases. As a result, 
the following conclusions can be made. 

1. There is a bias error in estimating rainfall from 
brightness temperature. 

2. This bias error is computable and is given by equation 
(3.4). The accuracy of computation depends on the place 
of truncation. 
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3. The bias eccor increases as the antenna resolution 
diecreases. 

Various methods for removing this bias are possible. The 
technique presented in the next section deals with estimating 
population variance i.e., (limiting value of variance as antenna 
resolution becomes higher and higher) from observed (or 
simulated) brightness temperatures* This technique has worked 
exceedingly well. 

3 .2 Error Due to Non-^Unigueness of and R Relat ion ship 

It is clear from Figure 2-1 that Tg and R relationship is not 
unique when estimating rainrate for observed brightness tem- 
perature, i.e., an observed Tg will result in two different 
values of R according to equation (2.4), Similarly, during 
simulation of ESMR data from GATE radar data two different rain 
rates may result in a single value for brightness temperature. 
For example, R=10 and R-10 9.5 mm/hr both map into Tg - 253, 6°K. 
Thus the simulated ESMR data inherently contains an error 
component which cannot be removed easily. Experimental results 
to indicate the magnitude of this error are presented in Section 
5. 
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4.0 UNBIASED ESTIMATION TECHNIQUE FOR MEAN imiNPALL 

This section deals with a mathematical formulation for the 
problem of finding an unbiased estimator of true mean rainfall 
based upon the brightness temperature measurements made by ESMR 
like instruments, the physical characteristics of ESMR involved 
in the derivation of Tq and R relationship, and adequate 
knowledge of the rainfall rate distribution. Also presented in 
this section are an estimation technique and a methodology for a 
practical implementation. 

4.1 Frequency Distribution of Rainfall Rate 

The frequency distribution analysis of GATE radar measured 
rainfall data was carried out for various fov's in the range 4km 
X 4km to 256km x 256km. For most fov's, this analysis indicated 
that rainfall rate, R, follows Gamma distribution. Figure 4-1 
illustrates an example for the case when fov is 8km x 8km. The 
rain rate Gamma density function, fj^(R, can be expressed 

as in equation (4.1), where and are the two parameters. 



It should be noted that both and are positive. 

The frequency analysis also showed that the rainrate exceeded 
20 mm/hr in less than 2 percent of the observations. As a result, 
for the sake of simplicity, one can truncate the Tg and R func- 
tional relationship of equation (2,4) at R in excess of 2 (j mm/hr 
so that equation (2.4) can be rewritten as; 
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Figure 4-1 

Rain Intensity Distribution for 8 KM Resolution 


Tg - F(R) - a - b exp (-cR) (4.2) 

Wherei a«271, b«107f and c«0.182 are constants. Knowing the 
general form of rain rate density function and the Tg and R 
functional relationship, one can find the density function of 
brightness temperatures, f^(TB r r PT)r follows (see 
Papoulis [1965]); 




^ ’ ^R * ?r2 

F’ (R) 


(4.3) 


Where P'(R) is the first derivative of F(R) in equation (4.2). 
It should be noted that f^i(TB , f ^t) also a Gamma density 
function with parameters Q’t and pijr both parameters being 
positive. Furthermore, from properties of Gamma distribution, 
the mean brightness temperature M.p (or T) and the variance ctq, 2 
can be expressed in terms of ^nd Pqi as given below. 


and. 


Similarly, 


T -^2 

15 = X 

“ cr.p2 

- 

"R 0-^2 


(4.4) 


and 



(4.5) 


where, (or R) is the mean rainfall and is the variance. 

The problem of estimating true mean rainfall can now be 
stated as follows: given an estimate of T and as obtained 


14 



from simulated ESMR data; the functional relationship between 
brightness temperature and rain rate; and thei form of rain rate 
distribution, find an estimate of mean rainfall, § (or Mj^), The 
mathematical details of this estimation technique are presented 
below. 


4.2 Mean Rainfall Estimation Technique 

The essence of the technique to be presented here is to 

A _ 

estimate R from the observed mean brightness temperature, T, and 
the variance, oij,2, of brightness temperature within a field of 

view. However, in order to facilitate easier understanding, Q'li 

and the parameters of rainfall rate distribution v/ould be 

estimated first. An estimate of mean rainfall can then be easily 

obtained from equations (4.5). 

The mean brightness temperature, T, can be expressed as 

00 

T “ J’'^B 

0 


Substituting for f«]i(TB , oip , from equation (4.3) and using 
the fact that dT = F'(R)dR, f becomes 



f R ( R J <*'r 3 Pr ) 

F' (R) 


F'(R) dR 


(4.6) 



(F(R) 


fR(R,Q'R,pR) dR 


Using the expression for P(R) from equation (4.2) and the 

» 

expression for fR(R, ^r , Pr) from equation (4.3) and simplify- 
ing, equation (4.6) results in the following. 
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T • a 


b 


Ml 

(pR+‘=)“‘‘ 


( 4 . 7 ) 


The intermediate eteps used in deriving equation (4,7) ate 
presented in Appendix A, 

The next step in the present estimation technique is to 
express in terms of o-j^ and The variance, can be 

expressed as E[T 3 ^J - T^, where E[Tg2j is the mean squared 
brightness temperature and can be expressed as in equation (4,8). 

r 


E [Tg ^ *J ^ 


(4.8) 


Once again, using the substitutions and simplifying in a manner 
similar to the ones used in deriving equation (4.7), we obtain: 


E [Tg^] = - 2ab 




+ b^ 


PR 


a 


R 


Pg+2. 


(4.9) 


The variance, criy^, can now be expressed as follows: 


0-^^ *= E [Tg^] - f ^ 



(4.10) 


Appendix A contains necessary steps used in deriving equations 
(4.9) and (4.10). 

Equations (4.7) and (4.10) express T and (r^2^ respectively, 
in terms of parameters and Pj^. Further simplifications, pre- 
sented in detail in Appendix A, result in the following; 
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and, 





I 






- Pr 




(pR+2cj 


0 


(4.11) 


(4.12) 


Where, and L 2 are functions of known quantities and are 
expressed as; 


In 


a - T 


(4,13) 


and. 


Jin 


OTn 


+ 1 


(4.14) 


I (a - T^) 

Note that T and are known since they can be computed from 
simulated ESf-IR data and a, b and c are constants used in equation 
(4,2) earlier. 

Equations (4,11) and (4.12) express and in terms of 
known quantities. However, unfortunately, these equations cannot 
be solved to give estimates of and pj^ in closed forms. As a 
result, one must use numerical methods to solve for and pj^ 
from these equations. An estimate of mean rainfall is then 
obtained from equations (4,5). 

It should be emphasized that the estimation technique 
presented in this subsection is valid for all fov's as long as 
the area of fov under consideration is at least equal to the 
spatial resolution of the antenna on the radiometer. 
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4,3 Analysis to Determine Spatial Cocrelation 

Spatial corcelation is an important parameter because it 
provides a measure of the characteristic rain cell size and, as a 
result, its determination will prove useful in the choice of 
adequate resolution for passive microwave radiometers. 

One way to formulate the problem is to consider the simulated 
ESMR data as a time series x(t), where x(t) represents a 
brightness temperature scene at time instant t. Each scene in 
this time series can be considered either as a spatial array of 
equal size cells or as a much longer linear array formed by 
taking consecutive rows or columns (from the spatial array) and 
placing them in a single row or column. For easier understanding 
of the mathematical details to follow, a linear array representa- 
tion is used. To remove the time variability, the time series 
x(t) is time-averaged to result in a single or master scene X, 
Note that this scene X may have high spatial variability. 

To study the spatial variability, the master scene, X, is 
scanned with different size antennas. Let represent the 
averaging distance of antenna under consideration, ESMR like 
radiometers, being integrating instruments, will provide indivi- 
dual brightness temperatures, TB(Da), which are distance apart 
from each other and can be expressed as in equation (4,15), 

T_(D^)=— / X(D) dD (4,15) 

° " Da i 
0 

Equation (4.15) emphasizes spatial variability by using X(D) to 
represent the master scene X, where D denotes distance. 
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It is cleat that the mean brightness temperature, Tg, of 
master scene is independent of since it represents the average 
of averages. However, the variance, is a more useful 
statistical parameter and is definitely dependent on It can 
be expressed as 







1 1 

(Da) E 



1 


1 

P3 

IB 

1 


Using equation (4.15) and simplifying! one obtains 


(4.16) 



(see Davenport and Root, p 69) 

Where Rjj(D') is spatial autocovariance function and is assumed to 
be exponential. Equation (4.17) can be further simplified to 
express o‘ii»^(D 3 ) as in equation (4.18). A detailed derivation of 
equation (4.18) is presented in Appendix B. 







^-1 + exp ( 



(4.18) 


Where, o-j^^ is the population variance or, alternatively, the 

D 

maximum value of Rj^fD) . And y = , Dq being the correlation 

distance, 
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The term correlation distance can be def ined in a number of 
ways. For example, it can be defined as that value of averaging 
distance, for which equals approximately 37 percent of 
its maximum value, R^Co). Alternatively, Dq may be defined as 
that value of for which equals approximately 73 percent of 

0 - 2 

X • 

The population variance and the correlation distance can both 
be estimated according to equation (4.18) on the basis of as few 
as two calculated values of o*(ii^(Da) computed for two different 
values of Dg. An unbiased estimate of true mean rainfall can then 
be obtained from and T 3 by using the p’\rameter estimation 
technique of subsection 4.2. This is discussed in more detail in 
the next subsection. 

4.4 Removal of Bias Due to Non-Line arity 

The estimation technique presented in subsection 4.2 works 
exceedingly well for all fov's and results in better estimation 
of mean rainfall than the one obtained by using transformation 
curve alone. However, it must be extended further in order to 
remove the bias in the estimate due to non-linearity of Tg and R 
relationship* 

As discussed in Section 3.0, the errors in the estimate of 
mean rainfall are caused by bias and are not random errors. As a 
result, statistical averaging over all fov's in the GATE size 
area is not very helpful. One approach to reduce or even 
eliminate this bias would be to use radiometers having higher 
resolution, i.e., capable of remotely sensing the brightness 
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tempo ratuces over smaller and smaller areas. As a matter of 
factf an ideal radiometer would be the one having the capability 
to measure brightness temperature at each point within the fov. 
However, improvements in radiometer resolution are limited by 
economical, engineering, and practical considerations since 
improvement in resolution translates into larger and larger 
antenna. 

A practical approach to removing the bias is as follows; 
using the radiometer of given resolution, find from observed data 
the mean brightness temperature, f, and the variance, oyer 
several larger fov's; plot crr^ versus the square root of the area 
of fov's ( k/A) and fit a multinomial function to obtain an 
analytical relationship between 0-^2 9^,3 ^nd extrapolate this 
curve to obtain an estimate of population variance t i.e., 
the value of cr^2 y;hen \/A approaches zero. The estimated popula- 
tion variance and T will then be used to estimate mean rainfall 
by the technique described in subsection 4.2. The resulting 

A 

estimate S will be unbiased since its expected value will be 
equal to the true mean rainfall. A detailed methodology for 
implementing this approach will be presented later in this 
section. 

An altei'/iate approach based upon the analysis of spatial 
correlation in subsection 4.3 has demonstrated good potential. 
In this approach, one needs to compute the variance, 
only two different averaging distances. Based upon these computa- 
tions, equation (4.18) is solved numerically, as described in 
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Appendix C, to give an estimate of population variance^ 
an estimate of correlation distance, Dq. The next step is to 
find an unbiased estimate of true mean rainfall. This is done by 
the same technique as in the previous approach. Some preliminary 
results obtained using alternate approach are quite encouraging 
and will be discussed in Section % 


4.5 Error Bounds on Estimate of Population Variance 

The estimate of population variance, as obtained by the 
technique described above will in general contain some error. 
This error is attributable, in la.'ge part, to the fact that the 
brightness temperatures can only be observed at discrete time 
intervals and not continuously. In otherwords, the sampling 
error can be used to determine error bounds on the estimate of 
population variance. 


A detailed analysis of sampling error has been performed by 
Laughlin and Gupta 11980] which resulted in the following expres- 
sion for sampling error, o-g^. 



Where, tg is the sampling period, tq is the temporal correlation 
interval, and T is the record length. The errors in estimating 
within each fov follow a normal distribution and, therefore, 
will be used to find error bounds to any desired confidence 
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level. For exampler within two standard deviations (20*^) the 
confidence level will be at 95%. As will be seen in Section 5,, 
these error bounds are used to find mean error in estimating mean 
rainfall as a function of antenna resolution. 

4r6 Methodology 

The methodology for obtaining an unbiased estimate of mean 
rainfall as described earlier can be implemented as follows: 

(a) Prom the simulated (or observed) datar compute the mean 

brightness temperature, T, and the variance, of 

brightness temperatures for several area fov's. Let A 
denote the area of an fov under consideration. Then 

^/A may be referred to as areal dimension or averaging 
distance. 

(b) Plot 0-^2 versus n/a and fit a multinomial function to 
obtain a functional relationship between the two. Note 
that f is same for all fovVs, since in each case the 
averaging is done over common area size, 

(c) Extrapolate the curve obtained in step (b) and obtain 

an estimate of population variance, This repre- 
sents the limiting value of <r^2 becomes diminish- 

ingly small. This will also provide an estimate of 
spatial correlation. 

(d) Use estimated population variance, and T to esti- 

A 

mate the mean rainfall, R, according to the estimation 
technique described in Section 4.2. R is an unbiased 
estimate of true mean rainfall. 
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(e) Use equ-xtion (4.19) to compute error bounds for desired 
confidence level. 

This methodology was implemented on the simulated ESMR data. 
The results obtained and a discussion of those results are 
presented in the next section. 


2A 


5.0 EXPERIMENTAL RESULTS AND DISCUSSION 


Earlier work of LaughXin and Gupta [1980] and other inves- 
tigators has shown that, from a sampling point of viewr samples 
taken every six hours are adequate. As a result, the simulated 
data from Phase I resulted in 72 brightness temperatures scenes 
and that from Phase II resulted in 64 brightness temperature 
scenes. From each time series, the mean brightness tempe|;ature, 
T, and variance, of brightness temperatures was computed for 

seven different size fov's: 4km x 4km, 8km x 8km, 16km x 16km, 

32km X 32km, 64km x 64km, 128km xI28km, and 256km x 256km. In 
each case the averaging was done over a common area of size 256km 
X 256km and, therefore, T remains same for all fov’s. However, 

does not and, in fact, as expected, it decreases as the area 
of fov increases. Computed values of T and are tabluated in 
Table 5-1. 

To determine the extent of errors due to non-linearity and 
non-uniqueness of Tg and R relationship, the mean rainfall rate 
for each Phase was computed using the Tg and R relationship of 
equation (2.4). The results are shown in Table 5-2, This non- 
parametric method was also used to find mean error in mean 
rainfall and the error variance which are shown in Figures 5-1 
through 5-4 for both Phases I and II. It is clear that both the 
mean error and the error variance increase with the size of the 
fov. Note that mean error is approximately 10 percent even for 
the smallest size fov. This error is due only to non-uniqueness 
of Tg and R relationship and non-linearity of this relationship 
does not affect it since no averaging was done for this fov. 
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TABLE 5-1 


MEAN BRIGHTNESS TEMPERATURE AND VARIANCES 
FOR PHASE I AND PHASE II 
FOR VARIOUS FIELD OF VIEWS 



Field Of View 

Phase 


4 km 

8km 

16 km 

32km 

64km 

128km 

r 

256km 

I 

Mean Brightness 
Temperature 
T 

168.6 

168.6 

168.6 

168.6 

168.6 

168.6 

168.6 

Variance 

267 

230 

190 

150 

105 

70 

30 

II 

Mean Brightness 
Temperature 
T 

167.4 

167.4 

167.4 

167.4 

167.4 

167.4 

167.4 

Variance 

0-^ 

198 

165 

126 

91 

55 

30 

16 
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TABLE 5-2 


TRUE AND ESTIMATED MEAN RAINFALL 
FOR PHASE I AND II 
OBTAINED BY NON-PARAMETRIC METHOD 


Phase 

True 

Mean 

Rainfall 

Estimated Mean Rainfall 

Field Of View 

4 km 

8km 

16km 

32km 

64km 

128km 

256 km 

I 

.468 

.443 

.379 

.339 

.309 

.286 

.269 

.251 

II 

.368 

.332 

.272 

.239 

.215 

.200 

.190 

.183 























Mean Error as % of 
True Mean Rainfall 



Figure 5-1 

Mean Error as % of True Mean vs. Areal Dimension 



100- PHASE II 

Mean Error as % of 
True Mean Rainfall 
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Figure 5-2 

Mean Error as % True Mean vs. Areal Dimension 
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This eccor component can not be easily corrected since it is 
inherently present in the simulated data itself. It should be 
pointed out that the true mean rainfall rate (0.486 mm/hr“^ for 
Phase I and 0.368 mm/hr“^ for Phase II) as obtained from radar 
measurements was assumed to be correct for the above method. 
However, as will be seen later, even the radar measurements were 
underestimates. 

The methodology described in subsection 4.6 was implemented 
using the values of T and f^om Table 5-1. The polynomial 
fits of sixth order were obtained for and were found to be 
very good approximations to the form in equation (4,18) which is 
expressed slightly differently in equation (5.1). 

+ - (- 1 + e 
A \ 




[a — ~ 

v/bere, has the same interpretation as in Section 4. The 
^ o o 

results are shown in Figures 5-5 and 5-6 for Phases I and II, 
respectively. The population variance, is estimated to be 

310 for Phase I and 23 0 for Phase II. The spatial correlation, 
is estimated to be 10km for Phase I and 8km for Phase II, 
The parameter estimation technique of subsection 4.2 resulted in 

A 

the unbiased estimate of mean rainfall rate, R, These results 
are tabulated in Table 5-3. The unbiased estimates are approx- 
imately 35 percent above the radar measurements which were made 
over 4km x 4km rain cells. This clearly demonstrates that even 
at 4km there exists a significant beam filling problem. At the 


32 









> 

to 

1 

\b 


£ 

C 

3 

§ 

g> 

Q 

UL 

1 


< 



I 


Variance at vs. Areal Dimension 


TABLE 5-3 


RESULTS FROM PARAMETER ESTIMATION TECHNIQUE 


Phase 

Population 

Variance 

Spatial 

Correlation 

Unbiased 
Estimate 
of Mean 
Rainfall 

I 

310 

10 

0.650 

II 

230 

8 

0.481 


TABLE 5-4 

RESULTS FROM ALTERNATE APPROACH 


Phase 

Population 
Var iance 

Spatial 

Correlation 

Unbiased 
Estimate 
of Mean 
Rainfall 

I 

308 

9 

0.641 

II 

226 

7.3 

0.462 
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same timer these are no potentiaXXy slgniticant gains to be made 
by using an antenna having a resolution of less than 1km or so. 

To test the performance of parameter estimation techniguer 

all computed values of were used to find corresponding 

estimates of mean rainfall at various antenna resolutions. The 

results are shown in Figures 5-7 through 5-10 which show the 

A , 

behavior of mean rainfall, R, corresponding variance, and 
standard deviation, This was done for both Phases I and II. 
The parameter estimation technique worked exceedingly well andf 
in addition, it was found that R, and all follow the 
functional form of equation 5.1 rather closely. 

Error bounds on the estimate of population variance, as 

obtained by evaluating equation (4,19) for six and twelve hour 
sampling rates are shown in Figures 5-11 and 5-12 for Phases I 
and II, respectively. These bounds are for 99% confidence level. 
Error bounds for one sigma confidence level were used to compute 
mean error (expressed as percentage of true mean rainfall) in 
estimating mean rainfall for various antenna resolutions. This 
was done for three measurement periods: 1 week, 2 weeks, and 1 

month. The results are displayed in Figures 5-13 and 5-14 for 
Phases I and II, respectively. As expected the mean error 
increases as the radiometer antenna resolution decreases, i.e., 
as the physical size the antenna is decreased. For example, an 
antenna having a resolution of 8km results in mean error of 
approximately 3 percent of monthly mean with two observations per 
day as compared to approximately 6 percent mean error for an 
antenna having a resolution of 32km. 
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Figure 5-9 

Variance a » and Standard Deviation a q vs. Areai Dimension 
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Figure 5-10 

Variance cr » and Standards Deviation Oq vs. Areal Dimension 



Figure 5-11 

Error Bounds on oi for Six and Twelve Hour Sampling Rates vs. Areal Dimension 



PHASE II 



Areal Dimension (KM) 

Figure 5-12 

Bounds on Ox Six and Twelve Hour Sampling Rates vs. Areal Dimension 
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Error Bound (Mean Error as % of Mean Rainfall) vs. Antenna Resolution 




PHASE II 


OF POOR QUALITY 
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Figure 5*14 

Error Bound (Mean Error as % of Mean Rainfaii) vs. Antenna Resolution 


The alternate approach developed in subsection 4.3 was tested 
and the preliminary results are fairly consistent with the 
results above. The correlation distance, was found to be 9km 
for Phase I and 7.3km for Phase II. These results are summarized 
in Table 5.4. This approach should be investigated further. 


6.0 SUMMARY 

The results of this work indicate that in order to obtain 
useful and unbiased estimates of monthly mean rainfall rate, an 
ESMR like instrument should be flown simultaneously on more than 
one satelite such that the effective re-visit time is approx- 
imately 6 hours. In addition, the resolution of the antenna 
should be consistent with the spatial correlation or correlation 
distance. The choice of antenna resolution of approximately 8 
kilometers for radiometers of future such as LAMMR* agrees quite 
closely with the spatial correlation of about 8 kilometers as 
found by this study. 

The bias due to non-uniform and incomplete filling of radio- 
meter beamwidth can be eliminated with reasonably high confidence 
by using the proposed technique. The unbiased estimate of mean 
rainfall rate over the GATE area was found to be approximately 35 
percent above the radar measurements. Also, the magnitude of 
error due to nonuniqueness of brightness temperature-rain rate 
relationship has been identified and found to be approximately 6 
percent of the unbiased estimate of true mean rainfall. 


* Large Aperture Muloi-channel Microwave Radiometer 
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APPENDIX A 
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APPENDIX A 

»lii?iww— IWMWHI tuiw— r ^ 

Derlvafcion of Equation (4.7) 

From equation (4,6), mean brightness temperature, T, is 
able as 


20 

f m I F(R) fj^(R,ag,3j^) dR 
0 


Substituting for F(R) and fg(R,aj^, , one gets 


T. 

J V 


, CR) V! ,<“R-^> ,, 


20 


m a 


$r"R 


■'(Or) 


R dR 


2-0 g OR -(S^+o)’' 

b / =; e 

'<»R> 


R 


(Otn-l) 


dR’ 


Now, using the fact that 

f j^(R 3|^) dR » 1 


/ 


we get. 


20 

I 




dR 


Hi 

Rr“R 


express- 


(A.l) 


(A. 2) 


(A. 3) 


A-1 
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Therefore, from equations (A, 2) and (A. 3), T, 


f 


a 


*^01r) (6r+c) 


“R 


becomes 


or 




m 


a 


gR°« 


This completes the derivation of equation (4.7) 


Derivation of Equations (4.9) and (4.10 ) 


From Equation 

(4.8), 


00 


MU. ' 1,1 M- 

E 

it/] 

SB 

/ 

f^ 

(Tr,c»,j,, 8^) dT 




J 

0 



Substituting T 

B=F(R) 

from equation (4.2) and for fm 

from equation 

(4.3), 

we 

get 



E [Tg^l - 

20 

f 

(a 

- b e“" 

cRj 

2 

J 




F’(R) 


a2 

20 

f 

Br“b 


($pR) (Cp-l) 

^ R HR 


C 

1 

^(«r) 


U£\ 


- 2ab 

20 

f 


e 

-(3p+c)R (oip-1) 

R R R ^ 



J 

‘<“r> 





20 

f 

3r“R 

T-* 

e 

(8p+2c)R 

dR 


F' (R) 


dR 


0 


^(Cr) 


A-2 


(A. 4) 


) 

dR 


ORIGINAL PAGE \S 

OF POOR QUALITY 


Now usiiig equation (A. 3) and simplifying, this becomes 


E [Tg^] » ~ 2ab ^ 


,ai 


3g+c 


+ b^ 


/ 3r 

I 3g+2c 


at 


(A. 5) 


Which is the same as equation (4.9). 


Now, variance, > can be expressed as ^ ■“ ^ 


s;2 


Substituting for the two terms on right and simplifying. 


a^ - 2ab 


Br \“r 


, 8 r+o 


+ b 


2 / PR 


O’, 


R 


3g+2cj 


a2 . b2 / 


2a 


R 


„ 2 k2 / 

0^ = b 


Br+o 


a, 


- 2ab 


^BR+cj 


B- \ ^“R 


Bg+2c, 


Br+c 


(A. 6) 


Which is the same as equation (4.10). 

Further simplifications to the above equations are presented 
below. From equation (A.4), 


or, 


3r\“R 


a-T 


b 


“r 


Br+c 


;=: £n| 


a-T 


or, 


a 


R 


[£n$g - £n(3g+c)] 


(A. 7) 
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Again, from equation (A. 4), b r 


?R 


in equation (A. 6) results in, 

2 




a T which when substituted 


tr,r 


(a-f ) 


4* 1 


(Pg+C) 


2a 


R 


an 

^ (a^+2c) 


or 


fn ( |--r + 1 


(a-T)2 


) - “r 


2Hn (3^+c) a. « Hn 


(3'j^+2c) 


0 *n 


\(a-T)Z ' *'""2 


■) 


Substituting for from equation (A. 7 ) and letting £n[ — I- g + 1 W L 
and simplifying, we get 

1^ 2jgn( 3j^+c) - fnfig - fn(3j^+2c) 

~ j2n(3j^+c) 


H H 


or 


or, 


(h ^H) - (4 •*• 4 ) *"Pr - 4«n(PR+2‘=) = 0 


£n 


(Pr-^4 


(1.2 + 2L^) n 


(L^ + L2> 




H 


0 


or, 


/ 2L-) 

(PrH 


(Lf + L2> 


= 1 


Pr 




or 


, . (L, + 2LJ (Lo + LJ, 

( - Pr ( 

which is the same as equation (4.12). 


0 


(A. 8) 
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APPENDIX B 

Derivation of Equation (4.t8) 


Equation (4.17) expresses »'rj, (D^) as in (B.l) 


dD' 


(B.l) 


2 “|D’I/D^ 

using, “ “*x ® 


O' 2 / ^ \ ^ 2 0*. 


/■(*-?) 


ID'I/D, 


dD' (B.2) 


a o 


The integrals involved in equation (B.2) can be easily evaluated 
to result in the following 


/ 


>?a 


D’ e 


dD’ = D^^ (l ^ e ^ °) - D„ D e ^ ° 
o \ / a 0 


and , 


D 


/ 


a -|D’|/D 

e ^ dD 


O jrv I 


°o( 


1 - e 


^D /D 
a' o 


D 

Substituting these in equation (B.2) and letting — 5. = y 

^o 


^(D„) = 


T 'a 


D 


D^U-e-Y) A. (l-e“Y) - D„ D^ e^^ 


D. 


a o 


2(r 


D 


1^0 - ^ (1 - e-y) 


B-1 




or; 


(T 2 




D D 

_e _ o (1 - e-y) 

°a ■'a 


or 


"■t^ <“a> “ 2<r/ 


J. 

L y 


•f ^ ( •" 1 4' 6 ^ ) 


(B.3) 


This completes the derivation. Further, in order to find the 

2 


limiting value of 0 ;^ (D^) as approaches 0, some further 
simplification of <B.3) is necessary as shown below. Also, 
L' Hospital's rule is used twice. 


him 


him 2 or 


*► 0 


y -► O 


Vjf (-1 + 


e-y> j 


Lim 2o* 


y o 


r i - e*"^ 
2y 


Lim 2 0 : 


y ^o 
2 

2 


•r [f] 


cr 


X 


That is, as the averaging distance becomes smaller and smaller, 


2 2 

tTm (D ) approaches , the population variance. 

■L cl X 


B-2 
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APPENDIX C 

Numerical Method to Estimate <r ^ and 

^ X o 


Let <r. ^ and 0 *-^ be the calculated values of <r 2 for averaging 

12 T 

distances D^ and D 2 , respectively. Further, let D 2 *= 20^^ for 
sake of simplicity. Then substituting these values in equation 
(4.18), one gets 


0-^^ « 2B 


and, 


■^ 2 ' = 


Dh I Dh 


D /D 


-1 + e 


-Di/Dq 


(C.l) 


2D, ^2D, 


-1 4* e 


•2D,/D^ 


Dividing equation (C.l) ty (C.2) and letting 


the following, 


or. 


0*, 


2 -I 


(C.2) 


k results in 


[Pj " Dq (- 


1 + e 


■D,/D 


°)1 


_2Di + (- 


-2D,/D 

1 + e 




= k 


(C.3) 


Let , 


Z = e 


-D^/Dq 


(C.4) 


Now, noting that (-1 + Z^) = (Z-1) (Z+1), equation (C.3) becomes 


4 [inZ - (Z - 1)1 _ ^ 
[ 2£nZ - (Z-1) (Z+1)] 


or, 


(4 - 2k) UnZ - 4Z + kZ^ + (4 - k) 0 


(C.5) 


C-1 


Equation (C.5) is solved numerically tor Z. The correlation 
distance, D^, is then obtained from equation (C.4) as 


D. 


o 


D. 


UnZ 


Substituting this in either equation 
will provide appropriate value of B. 


(C.l) or equation (C.3) 


